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Abstract 



Q^' The scattering theory of the integrable statistical models can be generalized to the 

case of systems with extended lines of defect. This is done by adding the reflection 
and transmission amplitudes for the interactions with the line of inhomegeneity 
^ ■ to the scattering amplitudes in the bulk. The factorization condition for the new 

5-( ' amplitudes gives rise to a set of Reflection- Transmission equations. The solutions 

of these equations in the case of diagonal S-matrix in the bulk are only those with 
5 = ±1. The choice S = —1 corresponds to the Ising model. We compute the exact 
expressions of the transmission and reflection amplitudes relative to the interaction 
of the Majorana fermion of the Ising model with the defect. These amplitudes 
present a weak-strong duality in the coupling constant, the self-dual points being 
the special values where the defect line acts as a reflecting surface. We also discuss 
the bosonic case S = 1 which presents instability properties and resonance states. 
Multi-defect systems which may give rise to a band structure are also considered. 
The exact expressions of correlation functions is obtained in terms of Form Factors 
of the bulk theory and matrix elements of the defect operator. 



1 Introduction 



In the past few years, methods and concepts of Quantum Field Theory (QFT) have been 
successfully applied to the analysis of homogeneous statistical models at and away from 
criticality. Besides the powerful techniques of Conformal Field Theories which enable us 
to disentangle the dynamics of the massless fluctuations of the critical points |1|, 0], the 
bootstrap approach - originally developed in particle physics to describe the scattering 
processes of strongly interacting particles - has proved to be one of the most efficient way 
to characterize the massive excitations of the statistical models away from criticality 0]. 
Exact solutions of massive QFT have been found along particular integrable trajectories 
of the Renormalization Group, characterized by the existence of an infinite number of 
conservation laws. Data easily accessible for the models associated to the integrable 
trajectories are provided by their exact factorizable ^-matrix and mass spectrum (see, for 
instance [3-8]). More importantly, the complete knowledge of the on-shell dynamics of 
those theories, combined together with general properties of analyticity and relativistic 
invariance, has allowed us to derive their off-shell behaviour, i.e. the computation of 
correlation functions, and to make contact with the original statistical models [9-15]. 

Homogeneous systems, however, are in many cases a mathematical idealization of the 
real physical samples which may present instead boundary effects and various types of 
inhomogeneity or defects. It is an interesting problem in statistical mechanics to estimate 
the influence of the inhomogeneities on the results obtained in pure cases and to develop 
the corresponding theory. With reference to systems with boundaries, they have been 
the subject of a wide investigation which has employed a large variety of techniques [18- 
22]. The bootstrap approach recently developed has brought new fight on the topic and 
has provided quite remarkable achievements in the understanding of QFT with boundary 
[23-29]. As we show in this paper, bootstrap methods are also extremely efficient to 
describe integrable statistical models with extended lines of defects. Before developing 
the bootstrap theory, it is worth to briefly discuss general aspects of statistical models 
with lines of defect in order to gain some insight to their properties [28-35]. 

One of the main reasons for considering extended lines of inhomogeneities is that only 
such kind of defects may affect the critical properties of the pure systems. Indeed, in 
the opposite case where there are only a finite number of localized inhomogeneities in 
the lattice, they would be eventually neutralized by iterating the Renormalization Group 
transformations so that the regime of the pure model will definitely take over. 

Scaling considerations are also useful to understand in simple terms the continuum 
version of the models with an infinite line of defect and to show that they may interpolate 
between a bulk or a boundary statistical behaviour. For sake of clarity, let us consider 
the simplest physical realization given by a system at temperature T in the bulk but 
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heated at a different temperature T along a line placed at the y axis. This system may 
be equally regarded as two semi-infinite copies of the model at temperature T coupled 
together through the energy density at the defect line. Its continuum properties are 
described by the euclidean action 

A = Ab + g J d^r5{x)e{r) , (1.1) 

where Ab is the action relative to the bulk and e(r) is the energy density with scaling 
dimension u. The scaling dimension of the coupling constant g = {T — T) is then given 
hy Ug = 1 — u. Consequently, all those systems with an irrelevant energy operator 
of scaling dimension u > 1 will exhibit the bulk critical behaviour near a defect line. 
On the contrary, those models which have a relevant energy operator with u < 1 will 
present a surface critical behaviour. The reason is that, in the former case the effective 
coupling constant may become arbitrary small and then the action reduces to that of the 
bulk theory, whereas, in the latter case it may take arbitrary large values suppressing 
all the fluctuations across the defect line between the two semi-infinite copies which will 
eventually decouple. 

An exception to the above pictures is given by the purely marginal case, i.e. u = 1 
which is realized in the Ising model. The interesting result obtained in the past by Bariev 
and McCoy and Perk ||3l| is that the model presents a non universal critical behaviour, 
with the critical indices of the magnetization operators continuously dependent on the 
parameter g of the action (|1 . 1| ) . The energy operator on the contrary remains a purely 
marginal operator for all values of the coupling constant g since its critical exponent u is 
fixed at the Ising value of 1 PUI ET, 



Let us now turn our attention to the boostrap theory of the integrable statistical mod- 
els with extended line of defect which was originally proposed in our previous publication 
[P7[] . In the continuum limit and away from criticality, such a theory can be formulated in 
terms of scattering processes of the massive excitations which take place either in the bulk 
or on the defect line. Hence, in addition to the bulk scattering amplitudes, we have to 
consider a new set of amplitudes relative to the interaction of the particles with the defect 
line. Because of the integrability condition, they only reduce to reflection and transmis- 
sion processes. It is then convenient to associate an extra operator D to the defect line 
and to formulate the dynamics in terms of algebraic relations which involve D and the 
operators A}^ of the massive excitations. The consistency of this algebra provides the 
unitary equations for the transmission and reflection amplitudes while the associativity 
condition gives rise to a set of cubic relations called Transmission-Reflection Equations. 
Although the general solution of these equations is still lacking, the Ising model with a 
line of defect can be identified as one particular solution of them. For this model, the 
particles in the bulk are given by massive Majorana fermions. The availability of an exact 
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resummation of the perturbative series in powers of the strength of the defect permits to 
test the complementary algebraic approach referred above. A similar approach is also used 
to discuss the case of free boson theory with a line of defect, which constitutes another 
solution of the Reflection- Transmission equations. The novelty of the model consists in 
the presence of resonance states and instability properties by varying the strength of the 
defect coupling. 

As for the integrable QFT in the bulk, we will show that also for the statistical models 
with a line of defect, the knowledge of the total set of scattering amplitudes permit the 
full reconstruction of the theory, i.e. the computation of the multipoint correlators. In 
particular, the aforementioned non-universal critical behaviour of the Ising model can be 
easily recovered by looking at the ultraviolet behaviour of the correlation functions. 

The layout of this paper is as follows. In the next section, we introduce the defect 
operator D and the relative algebra with the operators of the asymptotic particles A\. 
We initially derive the unitarity and crossing equations for the scattering amplitudes and 
then the Transmission-Reflection Equations which express their factorization properties. 
In Section 3, we discuss in detail one of the solutions of these equations, i.e. the Ising 
model with a line of defect. In Section 4 we consider the behaviour of the bosonic the- 
ory when coupled to a defect line. Geometrical configurations with richer structure of 
inhomegeneities are considered in Section 5. The computation of correlation functions of 
the model, based on the knowledge of the Form Factors of the theory in the bulk and on 
the matrix elements of the defect operator P, is presented in Section 6. Our conclusions 
are then summarized in Section 7. 



In the bulk, the theory of two-dimensional integrable statistical models with a finite cor- 
relation length can be elegantly formulated in terms of an ensemble of particle excitations 
in bootstrap interaction [3-8]. Although the bootstrap principle alone is in many cases 
sufficient to solve the dynamics, it is also useful to rely on a more conventional approach 
and to introduce an action that describes the interactions in the bulk 



In the following we assume we have completely solved the dynamics of the corresponding 
theory in the Minkoswki space and as a result, we know both the mass spectrum {m^} and 
the bulk scattering amplitudes S'^{[3ab) [|- The main features of the bulk scattering pro- 

^l^ab = Pa — (3b, where f3i is the rapidity variable of the particle It is related to the momenta by 
Pq"^ = rrii cosh/3i, p^*^ = rrii sinh/3i. 



2 Defect Algebra 




(2.1) 
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cesses can be briefly formulated as follows |Q . To each particle of the theory we associate 
an operator A^f]) and the set of all these operators satisfy the Faddev-Zamolodchikov 
algebra given by 

= S:t{f3ab) Al{(32) AiPi) ■ (2.2) 
The consistency of this algebra requires the validity of the unitarity equations 

s:tiP) si7i-P) = s'jr . (2.3) 



The commutation relations (p.2| ) should also be compatible with the requirements of 
algebraic associativity expressing the factorization of the scattering processes, i.e. they 
have to satisfy the Yang-Baxter equations 

SttiM Si:f,\{(3^,) Sil\l{(3,,) = ^ftt'(/?i3) Si%{(3,,) S^{[32s} • (2.4) 

The analytic continuation from the s-channel to the t-channel of the scattering process 
finally implies 

S'^m = ^|(^vr-/?) . (2.5) 

Eqs. ( |2.2| )-([275|), eventually supported by the bootstrap equations relative to the bound 
states 0, encode all the dynamics of the theory in the bulk. 

Let us now consider the presence of a defect line in the system, placed along the y-axis. 
The general action of the system is given in this case by 



A = Ab+ / d'r5ix)CD . (2.6) 




The new interaction, responsible for scattering processes which take place on the defect 
line, will generally spoil the original integrability of the theory: particles which hit the 
defect with sufficient energy may excite internal degrees of freedom of the defect (being 
eventually absorbed by it), or may give rise to production processes with multiparticle 
states propagating through the two semi-infinite systems placed on the two sides of the 
defect line. However, assuming that the additional interaction along the defect line is still 
compatible with the existence of an infinite number of conserved charges in involution, 
the dynamics drastically simplifies and consequently is suitable for an exact analysis, as 
we show in the sequel. 

By translation invariance along the ^/-direction (which we here identify with the time 



axis in the Minkowski space), for the theory described by the action (2^) we still have the 
conservation of the energy but not of the momentum. Therefore we may have scattering 
processes with an exchange of momentum on the defect line, compatible though with 
the conservation of the energy. If in addition to the energy other higher charges are 
also conserved, the scattering processes at the defect line must be completely elastic. In 
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particular, this means that a particle which hits the defect line with rapidity (3 can only 
proceed forward with the same rapidity or reverses its motion acquiring a rapidity of — /5. 
A further effect of the interaction with the defect line may be a change of the label of 
the particle inside its multiplet of degeneracy. The interactions of the particles \ (3;i > 
with the defect line will then be described in terms of the transmission and reflection 
amplitudes, denoted respectively by Tjj(/3) and Rij{P) (Fig. 1). 

The interaction of the particles at the line of inhomegeneity may be encoded in a set 
of algebraic relations analogous to those which describe the scattering processes in the 
bulk. In order to illustrate this, an additional operator D associated to the defect line 
should be introduced in the theory]^. This operator may be considered in relation to an 
additional particle state with zero rapidity in the entire time evolution of the system. Its 
commutation relations with the creation operators v4|(/5) associated to the asymptotic 
particles in the bulk are given by 



Al{(3)D = R,,{P)A]{-f3)D + T,^{P)DA]{(3) ; 
DAliP) = R,,{-(3)DA]{-P)+T,,{-(3)A]{(3)D 



(2.7) 



The first of these equations expresses the scattering of a particle that hits the defect 
coming from the semi-infinite system on the left hand side with rapidity f3. The second 
of ( |2.7| ) is obtained by an analytic continuation P — > — /5 of the scattering amplitudes of 
a particle that approaches the defect coming from the semi-infinite system on the right 
hand side. The consistency condition of this algebra requires the unitarity equations 

R.,i(3) R,,i-P) + T,,(/5) T,,i-/3) = ; ,^ . 

RmT,H-P)+TiiP)R^i-/3) =0. ^'^ 

Additional constraints emerge from the crossing relations 

R^Af^) = sfmRkiiiT^-P) ; .2 9) 



The first equation in (pl9| ) is obtained according to the argument proposed in p3|, p4| 
which exploits the quantization of the theory in the scheme where the time axis is placed 
along the x-axis. With reference to the second equation in ( |2.9| ), the transmission channel 
of the process shares the same properties of ordinary scattering in the bulk, the only 
difference being the occurrence of the particle D with zero rapidity. Thus, it is natural 
to assume that in the transmission channel the crossing symmetry is implemented in the 
usual way. We will assume the validity of eqs. ( p.9| ) and will check that they are actually 



^For simplicity, we discuss the case of a defect without internal degrees of freedom and therefore 
carrying no additional indices. The formulation of the more generale case is straightforward. 
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satisfied each time we will provide explicit solutions of the scattering theories with a line 
of defect. 

Usually the presence of an infinite number of integrals of motion implies not only 
the elasticity of all scattering processes but also their complete factorization, i.e. an 
n-particle scattering amplitude can be entirely expressed in terms of the elementary two- 
body interactions 0]. A crucial step for proving the factorization property of the total 
S'-matrix is to impose the associativity condition of the algebra ( p.7| ). In terms of physical 
process, this means that we prepare initially an asymptotic two-particle state consisting 
of I aI{Pi) Aj{P2) > with Pi > P2, and we let it scatter with the defect particle D with 
zero rapidity. The final output of the process should be independent from the temporal 
sequence of the elementary two-body interactions. Although what we have just described 
looks like an ordinary three-body process of the type that occurs in the bulk, there is 
however one distinguishing feature. In fact, in the three-body processes which take place 
in the bulk, given an initial state | Al{j3i)Aj{j32)Al{i33) > identified by a set of three 
ordered rapidities (3i > (32 > Ps, there is an unique final state given by the reverse ordering 
of the rapidities and possible exchange of the internal indices among the particles. On 
the contrary, for the scattering processes on the defect line we may have four possible 
final states, namely: (a) the state with both particles reflected by the defect line; (b) 
the state with both particles transmitted; (c and d) the states with one particle reflected 
whereas the other one transmitted. The final states may also differ from the initial one 
for the exchange of the internal indices and the above four possibilities give rise to a set 
of Reflection- Transmission (RT) equations shown in Fig. 2. 

The first of these (Fig. 2. a) coincides with the well-known boundary equations already 
analysed in |^, ^ 

Sg{Pi-P2)Rf,mS'^^{Pl+P2)RAf^l) = RaK{P2)Sit{P2+ Pl)Rf,{P2)Sl1{Pl- P2) ■ 

(2.10) 

The RT equations associated with the configurations of Figs. (2.b), (2.c) and (2.d) are 
given respectively by 

Si^iPi - P2) MPi) TUP2) = ^i1(/5i - P2) TUP2) T,i{Pi) ; 

StiPi - P2) TfbWi) Redm = Rcem SfJiPl + P2) TfbiPl) ; (2.11) 

S[^{Pl-P2)RfamS^^{Pl+P2)nb{Pl) = Tab{Pl)RM2) ■ 

Although a general solution of these equations is still lacking, it is easy to see that they 
become extremely restrictive once applied to QFT with a non-degenerate spectrum, i.e. 
those which have a diagonal S-matrix in the bulk. In fact, whereas eq. ( 2.10 ) and the first 
in ( p.llj ) are identically satisfied, the last two equations in ( 2.11 ) become in this case 

SabiPa + Pb) = SabiPb-Pa) , (2 12) 

SabiPa + Pb) SabiPa - Pb) = I , 
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Hence, from the first equation in ( p^.l2| ) we see that the S'-matrix in the bulk has to be 
a constant and from the second equation (or equivalently from the unitarity condition) 
this constant is fixed to be ±1. Thus we conclude that the only integrable QFT with 
diagonal S'-matrix in the bulk and factorizable scattering in the presence of the defect 
line are those associated to generalized-free theories. 

Obviously this restriction on the bulk S'-matrix does not apply when one considers 
purely reflective theories because they are ruled only by equations ( |2.10| ). Non-trivial 
solutions of these equations have been analysed for several models and they provide ex- 
plicit examples of QFT with boundary [23-29], some of them of relevant importance in 
statistical mechanics. 



3 Ising Model with a Line of Defect 

As we have seen at the end of the previous section, the validity of the Transmission- 
Reflection Equations in the case of non-degenerate mass spectrum selects S = —1 as a 
possible scattering matrix in the bulk. This solution can be identified as the scattering 
amplitude of the particle excitations of the Ising model, given by the massive Majorana 
fermions ^8] . The Lagrangian density of the continuum theory in the bulk is given by 



Cb = ^(x, t) (i7^a^ - m) ^(x, t) . (3.1) 

In the Majorana representation, given by 7° = (72, 7^ = —io^i, the fermionic field ^'(a;,^) 
is real, i.e. \E'^(a;, t) = \E'(x, t). The physical content of the model, as defined by the 
Lagrangian ( p.l| ), does not depend on the sign in front of the mass term since it can be 
altered by the transformation \E'— i>\E';\l'— \['of the fermionic field. As it is well known, 
the mass m is a linear measurement of the deviation of the temperature with respect to 
the critical one 

m = 2n{T - T^) , (3.2) 

and the symmetry m —m simply expresses the self-duality of the model. In the high 
temperature phase given by m > 0, the vacuum expectation value of the magnetization 
operator a vanishes, whereas, the corresponding quantity of the disorder operator yU is 
different from zero. Under the duality transformation, the role of order and disorder 
operators is reversed whereas for the energy operator e, given by e = i\Ef\Ef, we simply have 
a change in its sign. 

On a square lattice, the Ising model with a line of defect can be realized in two different 
ways (Fig. 3). The first is the chain geometry with bulk coupling constants J and modified 
couplings J parallel to the defect line. The second one is the ladder geometry, with the 
modified set of couplings placed in the perpendicular direction. Since the two geometric 
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realizations are related by Kramer- Wannier duality symmetry, from now on we can restrict 
our attention to one of them, say the chain geometry. In the continuum formulation, the 
defect line introduces the additional term^ 

Cd = -g6{x)'W{x,t)'^{x,t) (3.3) 

to the Lagrangian ( |3.1| ). The new interaction is purely marginal and therefore the beta- 
function associated to the coupling constant g is identically zero. The marginality of 
the interaction has as a consequence that the theory presents a non-universal ultraviolet 
behaviour in the magnetization sector, with the critical exponent of the magnetization 
operators which depend continuously on the parameter g, whereas the energy operator 



always keeps its original value of 1 of the Ising model [^O], q^, [33|| . 

In this section we are interested in determining the reflection R{f3) and transmission 
T{f3) amplitudes for the scattering of the fermion with the defect line, i.e. the S'-matrix 
elements between initial and final states u{pi) and u{pf) with pf = ±pi. To this aim, let 
us consider the perturbative series of the Green function of the fermion field based on 
the Feynman rules 

=^i2n)H\p-p')^^^^ . 9 . = - zg 2n 6{po - p'o) 
For the self-energy entering the exact propagator we have the following series of diagrams 



S V— = — + +»•»•>•> +••• 

where we have to integrate on the spatial component of the momentum running in the 
internal lines. The integral on the intermediate state is given by 

^= i-^g?^5{ko~po)J'i^,,^^ = -gHiko~po)^-^ • 
k 

In the above quantity we have discarded by parity the infinity related to the linear term 
in k. With this prescription, the geometric series for S is finite and can be expressed in 
a closed form as 

S(po) = 27vt6{po-Po) smx — : , (3.4) 

uj — im sm X 

where 



uo = ypo — "m? , sin x 



1 + f 



^The exact relationship between g and the lattice coupling constants will be established in section 6 
by comparing correlation functions computed in the lattice and in the continuum formulation. 
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We can now apply the usual LSZ reduction formulae, and for the transmission and reflec- 
tion amplitudes defined by 



out 



we have 



T{(3,g) 



cos X sinh /3 
sinh (3 — i sin x ' 

. sinx cosh/5 
1 

sinh f3 — i sin x 



(3.5) 



The transmission amplitude also contains the disconnected part relative to the free motion. 
Before commenting on the properties of these amplitudes, it is interesting to present 



an alternative derivation of (|3.5|) . This is obtained by implementing the algebra (|2.7|) on 
the creation operators of the fermion field. Let '^±{x, t) be the solutions of the free Dirac 
equation in the two intervals x > and x < 0, i.e. 



^(a;,t) = 9{x)'^+{x,t) + ^(-x)^_(x,t) 



(3.6) 



with the value at the origin given by ^1^(0,^) = | (\E'+(0,t) + \E'_(0,t)). The mode ex- 
pansion of the two components of the fields '$±{x,t) is expressed as 



V^«(x,t) 
^g)(x,t) 



2^ 



277 



iUe^ A(±)(/?) e-'^^it cosh 13-x sinh f^) ^ ^t^^(^) ^im{t cosh fS-x sinh 13) 

A(±)(/3)e-™(*™"''^-^'^^°^'^) + a;e-i^^)(/3)e^'"(*"°'^^-^' 



(3.7) 



with u = exp(27r/4), u = exp(— 'i7r/4). The operators A±{P) and A±(/3) satisfy the usual 
anti-commutation relations of a free fermion 



{A±{(3),Ai{e)} = 2n6{(3-e) 



although they are not all independent. They are related to each other by the conditions 
at a; = which arise from applying the eqs. of motion to ( p.6|) , i.e. 



(^f -^i'))(0,t) = f(^« + ^f3))(0,t) ; 
(^W_^«)(0,t) = f(^f + ^i^))(0,t) . 



(3.9) 



These equations are equivalent to the relationship between the modes 



M 



N 



Aiif3) 



(3.10) 
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where 



M 



N 



-i , g- E I , g H 

_^ o_ i. I 2 I 

cue 2 + -cue 2 we 2 uje 2 

9 g 

—uje2 — ^uoe 2 ujt 2 — |ci;e2 

§_ 2 S. . o_ i. 

—cue 2 — -one 2 — cje 2 + -cue 2 
9 9 



Hence, 




= M-^N ( •^-(-■'" ) = ( ) ( -^V"'' ) (3.11) 

AUP) ) [ T{,3,g) R(l3,g) ) [ aUP) j ' 

with R{P,g) and T{P,g) given in (|3.5|) . Note that, ahhough the boundary conditions 
( p.9|) are both hnear in g, there is however a feedback between the two components of 
the fermionic field. The final dependence from the coupling constant is then expressed in 
terms of trigonometric functions of the auxiliary angle x- 

Given the explicit expressions of the amplitudes ( |3.5|) , it is easy to check that they sat- 
isfy the unitarity and crossing equations (|2.8| ) and ( p.9|) . They present several interesting 
features. Firstly, by taking their sum and difference we obtain 



sinh 2 (/5 - ix) 

e'''^^TiP,g)-RiP,g) = 



2i5, - rr.f o ,.^ j^fo COsh i (/5 - 2x) 



cosh + ix) 

which can be considered as partial-wave phase shifts, with 5q and 5i crossed functions of 
each other. Secondly, notice that as functions of the coupling constant g, they satisfy a 
strong-weak duality given by 



-T{(3,g) , R[(3,- \ = R{(3,g) . (3.12) 





At the self-dual points g = ±2 the transmission amplitude vanishes and therefore the 
defect line behaves as a pure reflecting surface. From the unitarity equations (|2.8| ), the 
corresponding reflection amplitudes R{P, ±2) become pure phases, as can be explicitly 
seen by their equivalent expressions 

cosh I ± 2 sinh | 

m±2) = J^. . J • (3.13) 

cosh I =F I smh | 

They coincide with the reflection amplitudes of the Ising model with fixed and free bound- 
ary conditions respectively, as determined in [^. To establish directly the pure reflecting 
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properties of the defect line at the self-dual points, let us analyse more closely the decou- 
pling which occurs in the boundary conditions when g = ±2. For g = 2, the boundary 
conditions (B.9|) become 



(3.14) 



(^w_^(i))(o,t) = (v^f + v^i'))(o,t) , 

and taking their sum and difference, they can be written as 

(^(2)_^a))(o,t) = ; 

(^(2)+^a))(o,t) = . 

For g = —2, the original boundary conditions (|3.9|) are reduced instead to 

(^f +V^i!))(0,t) = ; 



(3.15) 



Equations ( p.l5| ) and ( |3.16| ) explicitly show that the two semi-infinite systems across 



the defect line are completely decoupled, and each of them can be treated as a QFT in 
the presence of pure reflecting surface whose role is to supply the appropriate boundary 
conditions [^|. At first sight, though, one may be surprised by the asymmetric form 



assumed by the eqs. (|3.15|) and (|3.16|) which treat differently the two fermionic fields 



'^±{x, t). However, this asymmetry has a physical origin. By means of the mode expansion 
(PT7|), the first equation in ( ^.15| ) and ( |3.16D can be used to determine directly the reflection 
amplitudes R{P, ±2). By the same token, using the second equation in (|3.15|) and ( p.l6| ) 



we find i?(— /?, ±2), instead. But, this is physically correct, the reason being that, in 
order to have a reflection of a particle described by \E'+(x, t) with the defect (boundary) 
line, this particle must approach the origin with positive rapidity (3. On the contrary, 
a reflection of a particle described by \I'_(x, t) with the defect (boundary) line is only 
realized for negative values of its rapidity. 

Further support of the identification of R{l3, ±2) with the reflection amplitudes of 
the Ising model with fixed and free boundary conditions comes from the analysis of the 
relationship between the lattice and the continuum formulations of the chain geometry, 
which will be established in Section 6. Anticipating the result, this is provided by the 
formula 

sinx = tanh2(J- J) . (3.17) 

Hence, the condition sin x = — 1 corresponds to a coupling constant J along the defect line 
infinitely larger (and positive) than the coupling constant J of the bulk. As a consequence, 
the spins along the defect line are frozen into a fixed boundary condition. On the other 
hand, the condition sin^ = 1 is obtained in the anti-ferromagnetic limit J — oo where 
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the spins along the defect hne are ahgned in antiparallel configurations. Since the nearby 
spins couple to a surface with vanishing magnetization, this situation corresponds to the 
free boundary conditions [p9|] . 

Let us now turn our attention to the analytic structure of the refiection and transmis- 
sion amplitudes. For negative values of the interaction with the defect line is attractive 
and consequently the theory presents a bound state with binding energy = mcosx- 
It is quite instructive to calculate the transmission and reflection amphtudes Tfe(/5) and 
Rb{P) relative to the scattering of the fermion with the excited state present on the defect 
line. The first thing to observe is that both amplitudes R{P) and T{j3) present a pole 
singularity at (3 = ix and P = ^(vr — x). The refiection amplitude -R(/3) has positive residue 
at both locations, given by isinx- On the other hand, T(/5) presents a positive residue 
with the same value as R{P) aX [3 = ix and a negative residue —i sin x at the other pole 
(3 = i{7c — x)- The problem of identifying which one of the two poles corresponds to the 
bound state is solved by selecting the singularity with positive residue in both amplitudes. 
This is the pole at /3 = ix- The relative binding energy is positive, as it should be. To 
recover the transmission and refiection amplitudes relative to the excited state, we have to 
impose the commutativity of the graphs shown in Fig. 4. Since the S'-matrix in the bulk 
is —1, the refiection amplitude Rb{P) coincides with the original one i.e. Rb{P) = R{P) 
whereas the transmission amplitude is given by Tb{l3) = —T{f3). If we again identify 
the singularity associated to a bound state as that pole with a positive residue in both 
channels, we see that for the defect bound state amplitudes the role of the two poles has 
been reversed! Namely, the pole which corresponds to the bound state in Rb{l3) and Tb{l3) 
is now located aX (3 = iij: — x) and is relative to the original ground state of the defect 



As a final remark of this section, the marginal nature of the defect interaction in the 
Ising model can be also inferred by looking at the high-energy limit of the amplitudes. 
For large values of (3 we have 

T(/3)~cosx, i?(/?)~2sinx . (3.18) 

Hence, except for the special values of the coupling constant g where one of the two 
quantities vanish, both amplitudes are always simultaneously present. Since this limit 
probes the short distance scales of the model, we see that its critical properties of bulk 



^Note that the presence of the transmission amphtude has been quite crucial in order to discriminate 
which one of the two poles with positive residue in the reflection channel corresponds to the bound state. 
In the pure reflecting situation, as for instance may be the case of the Ising model with a boundary mag- 



netic field considered in |23 , the occurrence of positive residue at both poles in the reflection amplitude 
and a misinterpretation of their role could in fact lead to a paradoxical hierarchy of bound states obtained 
by applying iteratively the boundary bootstrap equations. 
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and surface behaviour are simultaneously present]^. 



4 Bosonic Theory with a Line of Defect 



Another solution of the Reflection- Transmission equations is provided by the massive free 
bosonic theory with the S'-matrix in the bulk given by S = 1. As an example of a bosonic 
theory with a line of defect, we consider the model described by the Lagrangian 



(d^ifi) -m ip - g6{x)(p 



For the equation of motion we have 

□ + + g6{x)] ip = . 



(4.1) 



(4.2) 



As for the Ising model, one can obtain the reflection and transmission amplitudes by an 
exact resummation of the perturbative series in the coupling constant g. The calculations 
are analogous to the fermionic case, and rather than repeating them here, we prefer to 
exploit the algebraic approach directly. The solution of the equation of motion may be 
written as 

Lp{x,t) = 9{x)(f+{x,t) + 9{-x)ip^{x,t) , (4.3) 
where the mode decomposition of the two fields ^p±{x,t) is given by 

d(3 r 



^(±){x,t) 



27T 



j^^^^^p-j ^-im{t cosh l3-xsinh/3) _|_ ^^-j gim(i cosh /3-a; sinh /3) (4 4) 



The operators A±{f3) and ^±(/?) satisfy the usual commutation relations of a free massive 
boson 

A±(/3),4(e)] = 27r5i(3-e) . (4.5) 

The interaction along the defect however makes them not linearly independent. In fact, 
substituting eq. ( [4.3| ) into the equation of motion, the latter is equivalent to the boundary 
conditions 



d_ 

dx 



</p+(0,t)-vp„(0,t) 

iip4o,t)-ip4o,t)) 



; 

9 

4m 



(4.6) 



(¥.+ (0,t)+y.„(0,t)) 



^For an irrelevant interaction which leads to a bulk critical behaviour near the defect line, we expect 
in fact a vanishing of the reflection amplitude in the high-energy limit. On the contrary, for a relevant 
interaction, the system should show a purely surface critical behaviour characterized by the vanishing of 
the transmission amplitude. 
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which, in terms of the mode, can be written as 

R{f3,g) T{P,g) 

Ri(3,g) ' I ..^ I • i ■ ) 

The transmission and reflection amphtudes in the above formula are given by 

sinh f3 





sinh j3 + ig/ 4m ' 

ig/Am 
sinh (3 + ig/ Am 



T{(3,g) 
RiP,9) 

These amplitudes satisfy the unitarity and crossing equations ( p.8|) and ( p.9|) . It is easy 
to see that by substituting sinh (3 in ([4.8| ) with the linear momentum k, the two resulting 
amplitudes are the same as those obtained in one- dimensional quantum mechanics for 
the scattering in a 5-function potential (see, for instance, [0). However, due to the 
relativistic nature of the QFT, there is an important difference between the two cases, as 
shown by the analysis which follows on the pole structure of the amplitudes ( [4.8| ). 

For the 2iTi periodicity of the amplitudes, we can restrict our attention to the strip 
—ITT < /3 < in. Let us consider initially the case when g is a positive quantity. As long 
as g satisfies the condition < g < 4m, there are two poles on the negative imaginary 
axis relative to the unphysical sheet. By increasing the value of g they approach each 
other, and there is a critical value gc^ = 4m where they collide at position (3 = — i7r/2. 
Additional increment of the coupling constant causes the poles to move in the complex 
strip keeping their imaginary part equal to —in/ 2 but acquiring a real part (Fig. 5). In 
terms of QFT, this means that the bosonic theory with a coupling constant of the defect 
line larger than 4m presents two resonance states. As g grows, these poles move to 
infinity, and in the limit g oo, the defect line acts as a pure reflecting surface. Indeed, 
the transmission amplitude vanishes, whereas the reflection amplitude expresses the fixed 
boundary condition ip{0,t) = 0. 

Let us now analyse the case when is a negative quantity. In the range —4m < g < 0, 
the amplitudes present two poles placed on the positive imaginary axis relative to the 
physical sheet. The closest one to the origin can be interpreted as a defect bound state. 
By decreasing g, these two poles approach each other until they finally collide at (3 = in /2 
for the critical value gc2 = —4m. Further decrement of the coupling constant makes 
them move in the complex strip with an imaginary part equals to in/2 and with a real 
component which increases by decreasing g. However, these poles are now located in the 
physical strip and therefore the theory presents instability properties. The easiest way to 
explicitly illustrate this instability is to consider the analytic continuation /3 — > (*f ~ z^) 
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in R{(3). As discussed in Section 6, the resulting quantity R{(3), given by 

^^^^^ cosh P + g / Am ' ^ ^ 

can be interpreted as the amphtude relative to the emission of a pair of particles with 



momentum /3 and — /3 from the defect line placed along the x-axis |^. Then, for g < —4m, 
R{P) presents a pole for real values of (3 that induces a spontaneous emission of pairs of 
particles. The occurrence of such processes obviously spoils the stability of the theory. 
In light of the above results, we can summarize the discussion by saying that the 



QFT associated to the Lagrangian (4J.) makes sense only for values of g in the range 
—4m < g < oo. In a path integral approach to the problem, it is easy to see that 
there may be a competition in the Lagrangian (|4.1| ) between the genuine mass term and 
the defect interaction. Adopting the interpretation of the 5-function interaction as a 
suitable limit of a constant potential in the strip (— e, e) around the origin, when g is 
sufficiently positive in this interval, we may have an effective mass of the field if in this 
strip higher than the threshold mass m in the bulk. This produces the resonance poles 
in the transmission and reflection amplitudes. Viceversa, for negative values of g, the 
effective mass of the field ip in the tiny interval around the origin is smaller than the mass 
gap in the bulk and it decreases until it vanishes at g = —4m. After this value it becomes 
imaginary, giving rise to the instability property previously discussed. 

It is likewise interesting to understand the different behaviour of the bosonic and 
the fermionic theories in terms of the coupling constant. The reason is that the physical 
content of the fermionic model does not depend on the sign of the mass term, which enters 
linearly in the action. Therefore, by varying the coupling constant g, there is no a real 
competition with the genuine mass term in the action, so that the fermionic model cannot 
present instabilities or resonance states. In fact, crossing the critical values g = ±2, the 
poles simply interchange their positions, i.e. the weak coupling regime swaps with the 
strong coupling one. 

As a last comment on the bosonic theory analysed in this section, the defect interaction 
is associated to an irrelevant operator and therefore the defect line should be completely 
transparent in the ultraviolet limit. Indeed, taking the the high-energy limit /3 — oo of the 
amplitudes ([4. 81), the reflection amplitude vanishes whereas the transmission amplitude 
is identically equal to 1. 

5 Models with Multi-Defect Lines 

The solutions so far determined for the fermionic and bosonic theories in the presence of a 
single line of defect can be generalized and geometrical situations with a richer structure 
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of defect lines can be also included. In this section, we analyse the case of two parallel 
lines of defect, and then the quantization conditions induced by a periodic array of defects. 
Due to the different behaviour of the fermionic and bosonic theories, it is convenient to 
discuss them separately. 



5.1 Fermionic Theory 

Let us initially consider the Ising model with two parallel lines of defect, one placed at the 
origin along the y-axis with strength gi, the other shifted by a and with strength g2. In the 
fermionic formulation of the model, the field t) has a free motion in each of the three 
intervals /_ = (— oo,0), Iq = (0, a) and /+ = (a, +00) separated by the two defect lines. 
Therefore in each of the three intervals the field t) admits the usual decomposition in 
modes and the role of the defect lines is to provide the boundary conditions at the edges 
of the intervals. The first of them is at x = and is given by 



whereas for the second boundary condition at x = a we have 

(^f -4'^)(a,t) = f(^« + 4^))(a,t) ; 

(V^«-rf))(a,t) = f(^f + 4^))(a,t) . ^'-'^ 

In these equations the intervals are labelled by the subscript of the fields while their 
components by the upper indices. By using the notation Ri and Tj {i = 1, 2) for the 
reflection and the transmission amplitudes relative to the defect line with strength gi, it 
is easy to see that eliminating the intermediate modes relative to the interval /q, there is 
a linear relationship between the modes of the fields in the intervals /_ and /+ given by 



\ ^ 1 R{P,9i,92,a) T(p,gi,g2,a 
/ V T{/3,gi,g2,a) R{p,gi,g2,a 




(5.3) 



where 



T{P,gi,g2,a) 



R{P,gi,g2,a) 



l-riiP,a)R,{P)R2{P) ' 

R,{(3)+v{(3,a)R2mTm-R'M] 
l-r]{P,a)Ri{P)R2{P) 



(5.4) 



In the above expressions ?](/?, a) is a pure phase given by ri{P, a) = exp[— 2ima sinh P]. 
The above amplitudes satisfy the unitarity and crossing equations ( p. 81) and 



They describe the physical situation of a particle coming from the interval /_ with rapidity 
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(3 which hits the first defect line and, as result of this interaction, it can be either reflected 
or transmitted. When it is reflected, it appears as an asymptotic particle with rapidity —j3 
whereas when it is transmitted it approaches the next defect and can again be reflected or 
transmitted. As shown in Fig. 6, these two types of process may be repeated an arbitrary 
number of times at the two defect lines. 

Due to the existence of the fixed points g = ±2 of a single defect line, it is interesting 
to analyse some special limits of the expressions ( |5.4| ). To begin with, note that, at the 
values gi = ±2 where Ti = 0, the total transmission amplitude T{P, gi, g2,a) vanishes as 
well, whereas the reflection amplitude reduces to a pure phase given by R{P, ±2, g2, a) = 
±2). In this case, the first defect acts as a pure reflecting surface which therefore 
completely screens the presence of the second defect. The total transmission amplitude 
also vanishes when g2 = ±2. Concerning the reflection amplitude, it becomes a pure 
phase given by 

D/^ „ +9 - n(fH n\ R (R +9^ sinh (3{1 + r]-\(3, a) sin Xi) + i sin x(l - y-^P, a)) 
R{P,g„±2,a) - r^(/3,a)i?,(A±2) 3i^h/5(l + a) sinxi) - . sin x(l - a)) " 

(5.5) 

The total reflection process is now the result of an infinite sequence of elementary trans- 
mission and reflection scatterings at the first defect line combined with pure reflecting 
processes at the second defect line. Hence it is not surprising that the final expression 
depends on both R2{P, ±2) and the separation distance a. 

Except for the values of g when the defects behave as mirror surfaces, the possibil- 
ity for the fermion to go back and forth between the two defect lines produces typical 
resonance phenomena which are illustrated for instance by plotting the absolute value of 
T(/3, gi, g2, a). An example is shown in Fig. 7. 

Finally, by taking the limit a ^ 0, the two defect lines behave as a single one but with 
an effective strength g given by 

91 + 92 .p. „s 

This composition law of the defect strengths is similar to the addition of velocities in 
relativistic dynamics. The effective coupling constant g has as critical values g = ±2 
and reaches these limits when either gi or g2 are equal to ±2. This can be also seen by 
analysing the fixed points of the composition law defined by the iterative map 

9n + g 

9n+i = ^ 77 , (5.7) 

1 + gng/4: 

for some initial value g. 

The natural generalization of the situation with two defect lines is to consider a periodic 
array of defects all with equal strength g and separated by a distance a. The fermionic 
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field satisfies in this case the equation 



il^dfj_ — m — g ^ 5{x + no) 

n=—oo 



= 



and admits the decomposition 

oo 

m{x,t) = e{x -na)e{-x + {n + l)a)-^n{x,t) 



(5.9) 



with \I/„(x,t) solutions of the free Dirac equation. The dynamics of the model is entirely 
encoded into an infinite set of linear equations relative to the boundary conditions between 
the interval na and in — l)a, i.e. 



{A-€^){na,t) = f(V^i'i + V^«)(na,t) ; 
(^i'2i-^«)(na,t) = |(V;a + V;(2))(na,t) . 



(5.10) 



The simplest way to solve these equations is to employ a relativistic generalization of the 



Bloch theorem [40|, i.e to associate a wave vector k to the spinor field \1/ such that 



Equivalently, 



\E'„(raa,t) 



Aka 



*„_i((n- l)a,t) 



(5.11) 



(5.12) 



The wave vector k can always be confined to the first Brillouin zone —it /a < k < it /a. 
Plugging (|5.12|) into eqs. (|5.10|) , the resulting system is compatible provided that the 
equation 

sin (ma sinh j3) 



cos ka 



cosx 



cos (ma sinh f3) — sin 



(5.13) 



sinh P 

is valid. This equation gives rise to a band structure in the energy levels of the Majorana 
fermion of the Ising model, completely analogous to the periodic potentials considered 
in condensed matter physics. In fact, eq. ( |5.13|) can be satisfied for real k if and only 
if the right hand side of the equation is less than unity. Consequently, there will be 
allowed and forbidden regions of f3 and the corresponding spectrum of the energy, given 
hj E = m cosh (3, consists of a family of energy bands. A characteristic form of the 
spectrum is plotted in Fig. 8. For the pure reflecting values g = ±2, the above equation 
reduces to the quantization condition of the rapidity variable f3 



sinh (3 = ± tan(ma sinh (3) 



(5.14) 



which arises by considering the fermionic field defined in a strip of width a with fixed (+) 
or free (— ) boundary conditions at the edge of the interval. 
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5.2 Bosonic Theory 

The discussion of the bosonic theory largely follows the previous one and eqs. ( |5.4| ) is 
valid as it stands on the condition that we insert the bosonic amplitudes instead. Also in 
this case there are typical resonance phenomena produced by the trapping of the bosonic 
particle between the two defect lines. There is however a significant difference with respect 
to the fermionic case and this concerns the composition law relative to two defect lines 
with a separation a — 0. In this limit, the two defect lines behave as a single one with 
an effective strength g given by 

9 = 9i + 92 ■ (5.15) 

Due to the peculiar properties of the bosonic system discussed in Section 4, this compo- 
sition law implies that a system with two defect lines in the limit a —>■ may become 
unstable although each of the defect lines taken individually does not present any insta- 
bility property. Viceversa, one can obtain a well-defined bosonic system as a result of 
the limit a — > of a system which presents instability properties at one defect line and 
resonance states at the other. 

Taking the limit gi +oo, the first defect line becomes a pure reflecting surface and 
the total transmission amplitude vanishes. In this case the reflection amplitude reduces to 
R{l3,+oo,g2,a) = —1. The total transmission amplitude also vanishes when the second 
defect line acts as a pure reflecting surface. The corresponding reflection amplitude is a 
pure phase given by 

As in the fermionic case, the presence of an infinite periodic array of defect lines of 
strength g and separation a gives rise to a band structure described by a Kronig-Penney 
type equation 

7 f • 1. m , ^ sin(masinh/?) 

cos/ea = cosfmasmhpjH — — . (5.17) 

m smh p 

The pure reflective case g —>■ +oo gives rise to the quantization condition 

ma sinh /3 = nn , (ra = 0, ±1,...) (5.18) 

relative to the bosonic field in a strip of width a with fixed boundary conditions ip{0, t) = 
(p{a,t) = at the end points of the interval. 

6 Correlation Functions 

In the bulk, the scattering theory and the bootstrap approach -in addition to yield a clear 
understanding on the physical content of the continuum limit of the integrable statistical 



19 



models- have the added advantage of providing a powerful method for the computation of 
the correlation functions of the order parameters. Once the bulk S'-matrix is known, there 
are well-defined techniques for computing matrix elements of the local fields (j)i{x, t) on the 
set of asymptotic states < . . . /5„ | 0i(a;, t) | Pn+i, ■ ■ ■ Pm > and for reconstructing their 
correlation functions through spectral representation method based on the completeness 
relation of the asymptotic states. This program, known as Form Factor Approach, was 



originally proposed in |]^, |T0| and through this method many QFT have been recently 
solved [9-15]. In this section we will not go into the details of the Form Factor Approach 
in the bulk, which can be found in the aforementioned literature. We do however intend to 
prove that the spectral methods are also suitable for computing the correlation functions 
of the Ising model and the bosonic theory in the presence of a defect line. 

The easiest way to approach the problem is to use a formalism which takes full advan- 
tage of the solution of the theory in the bulk. To this aim, it is convenient to interchange 
the original role of the x and the t axes by the transformation x —it, t ix. The 
new space has a Minkowski structure with the defect line placed now at t = 0. In this 
new geometry, the space of the states is the same as in the bulk, and therefore, even in 
the presence of the defect line, the local operators 0j can be completely characterized by 
their known form factors. The presence of the defect line can be taken into account by 
defining an operator V placed at t = 0, acting on the bulk states. This operator plays 
the role of the S-matrix of the problem, and therefore, standard formulas of QFT allow 
the correlation functions to be expressed as |^| 



< 


T[(j)i{xi,ti) . . 


.V 


■ ■ ■ (j^niy'^ni t'n)\ 


> 


< 


V 


> 



In the above formula, are the fields in the Heisenberg representation, i.e. the 

representation where the time evolution is ruled by the exact Hamiltonian of the problem, 
including the defect interaction. On the other hand, are the field operators of 

the bulk theory and, as such, their time evolution operator is the bulk Hamiltonian^. The 
main advantage of eq. ( |6.1| ) is that, using the completeness relation of the bulk states, its 
right hand side can be entirely expressed in terms of the Form Factors of the bulk fields 
and the matrix elements of the operator V which are determined as follows. 

The defect operator T) encodes all information relative to the physical processes which 
take place at the defect line. To examine them, we have to initially realize that the first 
effect of the interchange of the x and the t axes consists in an analytic continuation of 



6, 



An equivalent way to look at eq. (6.1) is to consider a transfer matrix approach in the euclidean 



space. The transfer matrix may be written as T = exp[— i/^^] for all t but t = 0, where it is placed the 



defect line. Hence T) in (6.1) can be interpreted as the continuum limit of the transfer matrix operator 



which connects the states below and above the defect line. 
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the original rapidity (3 [i^ — the reason being that, to preserve the Minkowski 
structure in the new set of axes, we have to interchange correspondingly the momentum 
and the energy role. The rapidities are now measured as in Fig. 9. For convenience, it is 
useful to introduce the new transmission and reflection amplitudes, given by 

f = T - /?) , m = - . (6.2) 

They enter the expression of the simplest matrix elements of the operator P, given by 
Vi^i =< 13 \V \ e >, ©2,0 =< /9i,/52 I I > and ©0,2 =< Q \ V \ (3i,(32 >■ For the 
fermionic and the bosonic theory analysed in the previous sections, the first matrix element 
is easily computed by resumming the perturbative series with the defect interaction now 
localized at t — and the result is 

< l3\V\e>^ 27rf{l3)6{l3-e) . (6.3) 

By the same means, for the other two matrix elements, we have respectively 

</3i,/32 I P I 0> = 27rR{P,)S{P, + P2) , (6.4) 

and 

< I P I ^1, ^2 > = 27r^(^i) 5(61 + 62) . (6.5) 

Hence, T{l3) describes the process where a particle with rapidity j3 hits the defect line 
and is transmitted through it, keeping the same value of the rapidity (Fig. 9. a). On the 
contrary, R{P) may be interpreted as the amplitude for the creation or the annihilation of 
a pair of particles with equal and opposite rapidity (3 (Fig. 9.b). These three processes are 
compatible with the dynamics of the model because in a situation where the defect line 
is placed at t — 0, the processes are constrained by the conservation of the momentum 
but not of the energy. 

For the general matrix elements of the operator V, we can exploit the factorization 
property of the scattering theory and write down a set of recursive equations which involve 
the elementary two-body interactions considered above. For the bosonic case, the recursive 
equations are expressed by 

< (3i,...,(3i,...,(3m,P\'D\ei,...en>^ 
m 

= 27r ^ R{P)5{P + A) < . . . , A-i, A+i, ...Pm\'D\0i,...0n>+ (6.6) 

1=1 

+ 27rf2TW)s{f3-ej) <f3,,...Pm\T>\e,,...ej^i,ej+i,...,en> ; 

3=1 
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< Pi, ■ ■ ■ , Pm, \'D\9i,...9n,9> — 

n 

= 27r ^ R{e)5{e + 9^) <Pr,...,Pm\V\e^,..., 9i.i, 9^+,, ...,9n>+ (6.7) 



i=l 

m 



+27r T{9)5{9 - < A, . . . , . . . , /5„ | D | ^i, . . . , ^„ > . 

For the fermionic case, taking into account the anti-commutation relations of the fields, 
they can be written as 

<Pi,...,Pi,...,Pm,P\V\9u...9n>^ (6.8) 
= 27r 5^(-l)-+i- R{P)5{P + A) < A, . . . , A-i, A+i, . . . I P I ^1, . . . ^„ > + 

i=l 
n 

+2TiY.{-^yT{P)5{P-9j) <Pi,...Pm\V\9i,...9j^i,9j+i,...,9n> ; 

< /3i, . . . , Prn, I'D \ 9i, . . .9n,9 > — 

n 

- 27r Y.{-irR{9)5{9 + 9^) < P^, . . . , Pm\V \ 9,, . . . , 9,_,, 9^+,, . . . , > + (6.9) 

i=l 

m 

+27r ^(-l)'"-^r(^)5(^ - P,) <Pi,..., Pj-,, P,+i, ...,P^\V\9,,...,9r,> , 

These recursive equations can be graphically represented as in Fig. 10 and express the 
exact resummation of the perturbative series associated to the scattering matrix elements 
< m \ V \ n >. Since the particles are created or destroyed in couples, the non- vanishing 
matrix element < m \ V \ n > are only those with m — n — (mod 2). They are 
proportional to < | X> | > (which, for convenience, is set equal to 1) and the recur- 
sive equations permit to express all of them in terms of the elementary matrix elements 
^1,1)^2,0 a-nd 'Do,2 as previously determined. 

A useful method for solving the recursive equations is to introduce a generating func- 
tional of the matrix elements of T> by the formula 

g{ri, 7) = exp |/ dp ivi-PHP) + 7(/3)7(-/3)] + f{P)rj{Ph{P)j | . (6.10) 

Q depends on the two currents r){P) and j{P), which commute or anti-commute, depending 
on whether we are considering the bosonic theory or the fermionic one. The matrix 
elements of V are then given by 

d d d d 



< Pi,...,Pm:\T>\9,,...9n>^ (27r) — 



r)=7=0 

(6.11) 
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We are now in the position to compute correlation functions of local operators of the 
Ising model and the bosonic theory with a line of defect. Note that in computing the left 
hand side of eq. ( |6.1[ ) we should consider two different cases, namely: (a) the case where 
some of the operators are in the upper half-plane and the others are in the lower one, 
or (6) the case where the operators $j are all in one semi-plane, for example the upper 
one. In the former case, one has to use the general matrix elements < i \ V \ j >, and 
consequently both transmission and reflection amplitudes will enter the final expression 
of the correlation functions. In the latter case, on the contrary, the correlation functions 
will depend only on the reflection amplitudes R{f3) because, in virtue of the time ordering 
in eq. (|6.1| ), the defect operator V is in this case the last in the row and so, it acts directly 
on the vacuum state | >. Hence, the only matrix elements which enter the computation 
are Pj.o =< ^ | | >. Those describe the creation of the particle pairs and therefore 
only depend on R{P). 

In the remaining part of this section, using the form factors of the Ising model de- 



termined in H, 113] , and the matrix elements of the defect operator we compute some 



correlation functions of this model in the presence of the defect lineQ The simplest one is 
the one-point function of the energy operator e(x, t) which can be computed through the 
formula 

oo 

eo(t) = ^ < I e(a;,t) I n >< n I P I > . (6.12) 

n=0 

The energy operator couples the vacuum only to the two particle state, as can be easily 
checked by the fermionic representation of this operator, and for its matrix element we 
have 

< I e{x, t) I ^2 > = 27ri m sinh ^LZ_^ x (6.13) 
X exp [—mt (cosh [3i + cosh (32) + imx (sinh (3i + sinh (32)] , 

Hence the above sum (|6.12|) consists of only one term (Fig. 11) and using eq. (|6.4j) , it can 
be expressed as 

eo(t) = msinx H d(3^^^^ g-2mtcosh/3 _ ^g_^4) 

JO cosh p — sin X 

The one-point function does not depend on x, as it can be equivalently argued by transla- 
tion invariance along this axis. The above integral reduces to closed expressions in terms 
of Bessel functions when the defect line acts as pure reflecting surface. In the case of fixed 
boundary conditions, we have 

e{t) = -m [Ki{2mt) - Koi2mt)] , (6.15) 



'''All correlation functions will be computed in the euclidean space obtained by the analytic continuation 
t it. 
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whereas for free boundary conditions 

e{t) = m [Ki{2mt) + Ko{2mt)] , (6.16) 

In the general case, the one-point function interpolates between the two curves. The 
critical exponent of the energy operator in the presence of the defect line can be extracted 
by looking at the ultraviolet limit t — of its one-point function. For this limit we have 

/ N sin Y 

eo(t) ~ . (6.17) 

From this expression, we see that the defect line does not influence the critical exponent 
of the energy operator, which is the same as in the bulk, but rather enters the universal 
amplitude of the one-point function. For the pure reflecting case, the universal amplitudes 
coincide with those calculated in [Q. 

The relationship between the coupling constant in the continuum theory and in the 
discrete formulation can be extracted by comparing eq. (|6.17| ) with the analogous lattice 
computation, which reads 

, , tanh2(J — J) , , 

eo{t) ^ . (6.18) 

Hence, we have the following identification 

sinx = tanh2(J- J) . (6.19) 

In addition to the one-point function of the energy operator, it is also interesting to 
compute its two-point function. To simplify calculations, it is convenient to define the 
function 

2 J^oo cosh p — sin X 

Let us initially consider the situation where the energy density operators are on opposite 
sides of the defect line, i.e. ^2 > and ti < 0. The relevant expression in this case is 
given by0 

Giipi, p2) = E < I <P2) \ e(pi) I > . (6.21) 

id 

As before, the above series terminate. To explicitly evaluate it, in addition to the matrix 
elements ©2,0 and I^o,2, we also need the matrix element P2,2 given by 

+ f (/3i) T(/52)5(A - ^1) 6(02 - O2) + (6.22) 

- f(/?i)f(/32)5(/?l-^2)5(/?2-^l) 



^To simplify the notation, in the sequel we denote the couple of coordinate {xi,ti) simply by pi 
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With the notation t = t2 —ti and x = X2 — Xi, eq. ( |6.21| ) can be expressed as 



(- 



& 



d 



\dxdt 



dt 



+ 



(6.23) 



When the defect hne acts as a pure reflecting surface, all fluctuations across it are sup- 
pressed and this formula correctly reduces to the vacuum expectation values of the energy 
densities. 

Let us consider now the situation where the two energy operators are on the same 
side of the defect line, with t2 > ti > 0. For convenience, let us introduce the notation 



t = t2 — ti, t = t2 + ti, X = X2 — Xi and r = yxM-t^. The two point function can be 
written in this CclSG clS 



G2{pup2) = E < I e(P2) M >< ^ I e(pi) I J >< J I I? I > 



(6.24) 



There are only a finite number of non-vanishing matrix elements of the energy density 
and therefore the series truncates. It can be written as 



G2{pi,P2) = h + I2 + h 

where 

h = ^ r^i^^<0\eip2)\Pi,P2><Pi,P2\e{pi)\0><0\V\0> 

21 J~oo ZTC ZTT 

h - 
h - 

Ii coincides with the two-point function of the energy operator in the bulk. 



(6.25) 



hi 


+°° dl3i d/32 


< 1 




-00 2tc 271 


1 


/■+°° d(3i 


d[3i 


2! 2! 


J -00 271 


■ ~2^ 


1 


r+°° d(3i 


dPe 


2! 4! 


J -00 27[ 


■ 27r 



Ji = 



d 



d 



—KQimr) -\- —Koimr 



dx 



dt 



{Ko{mr)f 



The quantities which appear in I2 and I3 are the higher matrix elements of the energy 
density (which may be directly computed by the fermionic representation of this operator, 
e = i'^'^) and the matrix elements of the defect operator V, given by ( |6.11D . Considering 
that the computation of these quantities is lengthy but straightforward, we shall only 
present the final result 



I2 
h 



2 m sin x 



171^ sin^ X 



d 



d 



— iro(r)) {—Fix,t)]-Koir) 



^dx 



dx^ 



F{x,t) 



dx"^ 



\dxdt 



F{x,t) 



F{x,t) 



+ 



+ eo(ti)eo(t2) 
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Returning to eq. ( H.24 ), the two-point function can be cast in the form 
G2ipup2) 



eo(ti)eo(t2) + 
d 



- 2 




+ 





1 2 



dxdt 



F{x~t) 



Koir) + smx-g^F{x,t) 



1 2 



(6.26) 



It is now easy to verify that the expressions (|6.23|) and (|6.26|) coincide with those obtained 



in the lattice calculation 36 



As our last example, we discuss the one-point function of the magnetization operator 
cr(p) in the low temperature phase in the presence of the defect line. It can be calculated 
through the formula 



(To(t) = ^ < I a{p) \ nXn\V \ 0> . 

n=0 



(6.27) 



The magnetization operator couples the vacuum to all states with an even number of 
particles and its form factors are given by P, |l^ 



< I (t(0, 0) I a, . . . (32n > = (-0" n tanh : 



(6.28) 



Since the matrix elements of V in ( |6.27|) are different from zero only for pairs of particles 
of opposite momentum, we are lead to consider the matrix elements of the magnetization 
operator given by < | cr(0) | — . . . , —Pn,Pn >■ They can be conveniently written 



as 



< I a(0, 0) I -pu A, ... - pn, /5n > = ^" [^H tanh AJ x det W{p^, pj) , (6.29) 

where W{Pi,Pj) is the n x n matrix given by 

2 J cosh Pi cosh Pj \ 

(6.30) 



WiP.,Pj) 



cosh Pi + cosh Pj 



Hence, the one-point function is the sum of an infinite number of terms shown in Fig. 12 
and it can be expressed Fredholm determinant 

^o{t) = E ^ / dp,...dPn{ll^ tanh A R{P,) e-^-^^^^A detiy(A, Pj) = 

n=0^'-^-°° \k=0 J 

= Det{l + zW) . (6.31) 



The explicit form of the kernel is given by 

E{pi,mt,x)E{pj,mt,x) 



>v(A,A-,x) 



cosh Pi + cosh Pj 



(6.32) 
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where 



E{/3,mt,x) = sinh /?e-™*'=°^^^ (cosh /3- sin x)"^/^ , z = . (6.33) 

ZTT 

In terms of the eigenvalues of the integral operator and their multiplicity, o"o(t) can be 
also expressed as 

oo 

Mt) = Ui^ + ^^^r (6.34) 

i=l 

As far as mt is finite, the kernel is square integrable and therefore all results valid for 



bounded symmetric integral operators apply (see, for instance [0). In particular, for 
mt oo, crQ{t) falls off exponentially to the bulk vacuum expectation value. However, 
when mt 0, the integral operator becomes unbounded. The multiplicity of the eigen- 
values grows logarithmically as a ~ Mn ^ whereas the eigenvalues become dense in the 
interval (0, oo) according to the distribution 

A(p) = . (6.35) 
cosh Tip 

Hence, for the critical exponent of the magnetization operator, defined by 

Mt) - , (6.36) 



we have 

Xa{x) = - 



1 r°° [ 2tiz \ 11 

- / dp In IH — = -- + — — arccos^(- sinx) • (6.37) 

TT Jo \ COShp j 8 ZTT^ 



This expression agrees with the lattice calculations [^, ^ and since it depends on the 
coupling constant, it explicitly shows the non- universality of the model. 



7 Conclusion 

The main purpose of this paper was to prove that the bootstrap approach can be success- 
fully extended to integrable models with linear inhomogeneities and that the computation 
of the correlation functions for those systems can be achieved by means of a suitable gen- 
eralization of the Form Factor techniques. We have analysed the general situation in 
which translation invariance is broken by the presence of defect lines allowing reflection 
and transmission processes. While at the moment it is still an open problem to see 
whether there are other solutions of the Transmission-Reflection equations in addition to 
the fermionic and the bosonic theories analysed in the text, it is worth to stress that the 
method for the computation of the correlation functions exposed in Sec. 6 is expected to 
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work without limitation in the pure reflecting case. This corresponds to the boundary 
field theories which have recently received a lot of attention in view of their potential 
application to a wide class of physical situations. 

Acknowledgements. We are grateful to R. lengo and M. Nolasco for useful discussions. 
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Figure Captions 



Figure 1 . Reflection and Transmission Amplitudes. 

Figure 2 . Reflection- Transmission Equations. 

Figure 3 . Chain and ladder geometry of a defect line lattice. 

Figure 4 . Bootstrap Equations of the defect bound state in the reflection (a and b) and 
in the transmission channel (c and d). 

Figure 5 . Pole structure of the bosonic amphtudes for positive values of the coupling 
constant (empty circles) and for negative ones (flUed circles). 

Figure 6 . Scattering processes at the two defect hues. 

Figure 7 . Resonances in the transmission channel of the scattering on the two defect 
lines. 

Figure 8 . Band structure of the energy levels of the Major ana fermion with an inflnite 
array of defects. 

Figure 9 . Defect line at t = 0. Transmission channel (a) and processes of creation 
(annihilation) of a pair of particles {b). 

Figure 10 . Recursive equations of the matrix elements of the operator T>. 

Figure 11 . One point function of the energy operator of the Ising model in the presence 
of the defect hne. 

Figure 12 . One point function of the magnetization operator of the Ising model in the 
presence of the defect hne. 
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